Is the PAMELA anomaly caused by the supernova explosions 

near the Earth? 
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We show that the anomaly of the positron fraction observed by the PAMELA experiment can 
be attributed to recent supernova explosion(s) in a dense gas cloud (DC) near the Earth. Protons 
are accelerated around the supernova remnant (SNR) . Electrons and positrons are created through 
hadronic interactions inside the DC. Their spectrum is harder than that of the background because 
the SNR spends much time in a radiative phase. Our scenario predicts that the anti-proton flux 
dominates that of the background for > 100 GeV. We compare the results with observations (Fermi, 
HESS, PPB-BETS, and ATIC). 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 



Observations with the PAMELA ^ and ATIC/PPB- 
BETS [2] experiments have imphed or shown anomaly 
in the positron fraction and a sharp feature in the CR 
electron spectrum, although Fermi satellite recently re- 
ported that the feature in the electron spectrum is much 
smaller or even absent depending on the uncertainty of 
the background cosmic-ray electron and positron spec- 
trum 0. The PAMELA (and ATIC/PPB-BETS) results 
indicate that there should be an unknown source of elec- 
trons and positrons. One attractive idea of the source is 
the annihilation or decay of weakly interacting dark mat- 



ter particles 0, 
such as pulsars 



1. The other is astrophysical scenarios 



Among the astrophysical sources, supernova remnants 
(SNRs) are one of the promising candidates for the elec- 
tron and positron source. Their synchrotron emission 
reveals the existence of high energy electrons in SNRs 
Q , and their 7-ray emission would indicate that of high 
energy protons [9] , although there has been no direct ev- 
idence of proton acceleration in SNRs yet. 

In this paper, we show that the observed electron and 
positron excesses can naturally be explained by recent 
supernova explosions near the Earth. In fact the low- 
density local bubble (LB), in which the solar system lo- 
cates, and the adjacent bubble (Loop I) with sizes of 
~100 pc are often thought to be created through ^20-40 
supernova explosions that started ~ 10'' yrs ago [lOj . 
Since those explosions occurred close to us (~ 100- 
200 pc) , they could have affected the environment around 
the solar system more significantly than the supernova 
explosions in the spiral arms '--^kpc away pT| . 



II. SUPERNOVA(E) IN A DENSE CLOUD 

We consider supernova explosions that happened ~ 
10^~^ yrs ago in a dense gas cloud (DC) ~100-200 pc 
away from the Earth. The DC does not necessarily cor- 
respond to the host DCs of the supernovae that created 
the LB or Loop I. Massive stars that explode as super- 
novae tend to be born in giant DCs |12l |. They may 
explode near or inside the host DCs, because their du- 
rations of life are short. In this paper, we assume that 
although a DC had not been destroyed by ultra-violet ra- 
diation from massive, short-lived OB stars, it is ionized 
and the temperature is ~ 10"* K when those stars explode 
[l^ . In such environments, the shock front of an SNR 
accelerates protons, which create electrons and positrons 
through hadronic interactions in the surrounding DC. 

The evolution of a SNR and the particle acceleration 
around it is modeled following a previous study (14i | . The 
shock velocity Vg is written as a function of the SNR age 
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where Vi — Vi^glO^ cm s ^ is the initial velocity, and 
SsN = Ec^ilO^^ erg is the initial energy of the ejecta. 
When we consider multiple supernova explosions, we 
simply assume that they explode almost simultaneously 
(within a time-scale of the destruction of the DC) and 
give a larger Esn than that of a single supernova. At 
tago < ti = 'l5{E5i/n2y^^v~^^^ yr, the SNR is in the free 

expansion phase, at tago > t2 — 3.5 x lO'^E'g^^^rtj ^^^^ yr, 
it is in the radiative phase, and at ti < tagc < ^2, it is in 
the Sedov phase. Here no = ^2 10^ cm^'^ is the proton 
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density of the DC. 

We assume that the energy spectrum of the accelerated 
protons is 



Np{E) cx E-^expi-E/En 



(2) 



The index is represented by s = (r + 2)/(r — 1), where r 
is the compression ratio of the shock [l5| . which is given 
by the Rankine-Hugoniot relation for ti < iagc < a-nd 
r ^ \/2vs/vAu for iagc > ^2, where vau is the upstream 
Alfven velocity 14]. The maximum energy is determined 
by the age of SNR: 
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where Vs,d, — t's/lO^cm s ^ ^ h 1) is the factor deter- 
mined by the shock angle and the gyro-factor, and i?d 
is the downstream magnetic field, which is expressed by 
= fSoCj where Bdc is the magnetic field in the DC 
[l^ . Since we consider an ionized DC, we ignore the ion- 
neutral wave damping at the upstream of a shock (for 
particle acceleration in a radiative phase, see [IBl)- We 
equate the minimum energy of the protons with the pro- 
ton rest-mass energy. 

We assume that supernovae explode at the center of a 
DC for the sake of simplicity. We take Vi^g = 1, /i = 1, 
and i?DC = 10 fJ-G from now on. 

We consider the case where a few supernovae (or one 
hypernova) explode almost simultaneously. Thus we as- 
sume that E^i = 3, from which we obtain ti = 65 yr and 
t2 = 4.5 X 10^ yr, if the proton density inside the DC is 
no = 100cm"3. At iagc = 8.0 x lO'' yr, the Mach num- 
ber and radius of the SNR are Ala = 7 and i? = 26 pc, 
respectively. At this point, the SNR is in the radiative 
phase. We assume that particle acceleration around a 
SNR ceases for Ma < 7 (see e.g. [l7|), although the 
exact value is not clear. We refer to the age at which 

Ala — 7 as ^agc — ^acc ■ 

At iago = tacc, we obtain s = 1.0 and i?max,p ~ 
120 TeV. Although particle acceleration stops at this 
time, the SNR continues to expand. The radius is 
Rs — 50 pc at iagc = 5 X 10^ yr. Since it is compa- 
rable to the size of a giant DC [3 and Esn is larger 
than the binding energy of a DC, the cloud would be 
destroyed around this time. Until destroyed, the DC is 
illuminated by the accelerated protons from the inside 
with the spectrum of Eq. ^ given at iagc ^ ^acc- The 
duration of the exposure (tpp) could be approximated by 
the time elapsing from the explosion of the supernovae 
to the destruction of the DC because of the short time- 
scale of iacc- The protons create e~, C*", p and p through 
pp collisions in the DC for the period of ~ tpp. Based 
on this scenario, we calculate their spectrum observed 
at the Earth, and search parameters that are consistent 
with observations. Since the model of the SNR evolution 
may have some uncertainties, we do not stick to the exact 
values of parameters in the above discussion. 



After the destruction of the DC, the created e~, e"*", p 
and p propagate through diffusion processes and reach to 
the Earth. Since we assume that the DC has already been 
destroyed until the current epoch, there is a difference in 
the last arrival time between massless neutral particles 
such as photons and neutrinos, and charged particles. 
We would not detect any photon and neutrino signals 
from the DC ~ 10^~^ yrs after the destruction. 

We have calculated spectra of those daughter particles 
of the pp collisions by performing the PYTHIA Monte- 
Carlo event generator [l^. Then we solve the diffusion 
equation of the charged particle "z" (i runs e~, e^, p and 
P), 



me,)h] + g(£,) (4) 



where fi(t,x,ei) is the distribution function of i and 
£j — Ei/GeV with Ei being the energy of i. The flux 
^i{t,x,ei) is given by <i>i = (c/47r)/i. For electrons and 
positrons, according to a model 08-005 in ^205, we take 



X(£e) = i^o(l+ee/3GeV)^ 



(5) 



with Ko ^ 2 X lO^^cm^s"^ and S = 0.6. The coohng 
rate through the synchrotron emission and the inverse 
Compton scattering is given by 

B{ee) - 10-^%-^el[0.2{Bdis/3fiGf + 0.9]/l.l (6) 

[2ll |. where Sdiff is the magnetic field outside the DC. 
These parameters approximately correspond to so-called 
"med model" of the particle propagation . 

If we assume that the timescale of the production is 
shorter than the diffusion timescale ~ cP/{Kc), and the 
daughter-particle source is spatially localized well, we can 
also use the analytical solution in ^2^]. When the shape 
of the source spectrum is a power-law with an index a, 
Q = Qo£~°'^{x)S{t), the solution is represented by 
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with £cut = (^^diff) ^ and the diffusion length 



1 - (1 - Ee/ecut)!-*' 
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Here we took an effective distance to the source d, which 
is calculated by spatially averaging the distance to the 
volume element of the source, and a ~ s in the current 
calculation. The time after the input of the daughter 
particles is defined by idiff- Approximately QqS'" is rep- 
resented by Vstppd'^rii/ {dtdEi) with Vg the volume of 
the source and 



dtdE, 
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Here the differential cross section of tlie "j"-niode for the 
production of i particle is daj{Ep, Ei)/dEi with its mul- 
tiplicity gj = gj {Ep ,Ei). We also consider the neutron 
decay for the electron/positron production, which is not 
included in the original version of PYTHIA. The initial 
proton spectrum Np{Ep) is normalized to be 

Vs J dEpNp{Ep) = Stot,p , (10) 

which is given later. 

For the local propagation of p and p, the cooling term 
is negligible. In addition, we may omit the annihilation 
of antiprotons through scattering off the background pro- 
tons because the mean free path for the annihilation re- 
action is large. We can also omit the convection by the 
convective wind. Thus the analytical solution for the 
proton and antiproton can be also given by Eq. ([7]) with 
taking a limit of ep/ccut — *■ 0. 

III. RESULTS 

Fig. [1] shows the positron fraction and the e^-e+ flux 
observed at the Earth. We assume that a spherical DC 
with a radius of i?Dc — 40 pc and proton density of 
no — 100 cm^"^ is uniformly illuminated by protons with 
the spectrum given by Eq. (U) with s — lA and -Emax = 
100 TeV for the duration of tpp ^ 5 x 10^ yr. The total 
energy of the primary cosmic-ray protons is taken to be 
£"101, p = 3 X 10^° erg. The distance to the front side of the 
DC is d = 200 pc. We assume that the diffusion time of 
and e+ is idiff = 5 x 10^ yr, and Bdiff = 3/iG. We refer 
to this parameter set as Case A. Observational results are 
also shown. This result is consistent with observations for 
the positron fraction, and the e~-e'^ flux obtained with 
PPB-BETS and ATIC2. 

The feature of the positron spectrum can intuitively 
be understood as follows. The spectrum of the secondary 
e~ and e+ from the DC is harder than that of the back- 
ground one, because the primary protons are mainly pro- 
duced in the radiative phase of the SNR with the com- 
pression ratio r much larger than 4. The cutoff scale is 
given by E,^t = (^tdiff)"' GeV - 600 GeV/(idiff/5 x 
lO^yrs). The sharp rise around the cutoff energy is made 
by and e+ accumulated through cooling. That is be- 
cause we are considering the case of s < 2, which means 
that many of the cooled electrons sizably contribute to 
modifying the spectral shape around the cutoff energy. 
However, the actual spectrum would be broader with a 
width of Ai?cut /-E-cut ~ tpp/tdis ~ 1, because the injec- 
tion continues for a finite duration tpp 

Fig- m is another example of the positron fraction and 
the e~-e+ flux. In this case (Case B), we adopt i?DC = 
40 pc, no = 50 cm-^, s = 1.75 , E,^^^ = 100 TeV, tpp = 
2 X 10^ yr, Etoup = 3 x 10^" erg, d = 200 pc, tarn = 
2 X 10^ yr, and B^m — 3 /iG. The results are consistent 
with observations for the positron fraction and the e~-e~^ 
flux obtained with Fermi. Since we have not performed 
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FIG. 1: (a) Positron fraction for Case A (solid line), which 
includes the electrons and positrons coming from the DC and 
background electrons (dotted line). Filled circles correspond 
to the PAMELA data (b) Electron and positron flux for 
Case A (solid line). The flux of the electrons and positrons 
created in the DC is shown by the dashed line, and that of the 
background ones is shown by the dotted line. Observational 
resuhs (Fermi, HESS, BETS, PPB-BETS, and ATIC2 [1 i, 
^^) are also presented. The shadow shows the uncertainty of 
the HESS data. 
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FIG. 2: Same as Fig. [1] but for Case B. 



a complete parameter search. Cases A and B should be 
regarded as examples among possibilities. 

Fig. [3] shows the antiproton fraction at the Earth. In 
our scenario, the flux of anti-protons dominates that of 
the background for > 100 GeV, which can be used to dis- 
criminate our hadronic scenario from other leptonic ones. 
Antiproton fraction rises at higher energies than that of 
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FIG. 3: Antiproton fraction for Cases A and B (solid lines). 
The dotted line shows the fraction calculated only from the 
background flux. The fitting formula of the background pro- 
ton and antiproton are taken by [2^. We adopted data points 
taken with PAMELA [H. 



the positron fraction (Figs. [T] and H]), because the spec- 
tral index of the background antiproton is harder than 
that of the background positron. So far the PAMELA 
antiproton data could not have distinguished one model 
from another yet (Fig. [3]). 

We note that old SNRs in DCs could be observed as 
unidentified TeV sources [3, [13, . Broad-band spectra 
of the old SNRs could constrain the parameters of our 
scenario. It is likely that the DC has made a star cluster, 



which might be identified near the Earth. 

Recent observations have shown possible decline of 
B/C ratio toward higher energies (> 100 GeV), which 
seems to be inconsistent with the present scenario be- 
cause SNRs also accelerate metals |29l . Isoj . In our sce- 
nario, however, the source of electrons and positrons is 
the SNR in a specific dense cloud, and it is very localized. 
Thus it is subjected to the variation of metal abundances 
in the Galaxy. In fact the abundance of Carbon in Ca- 
rina Nebula, which is a star forming region in the Galaxy, 
is < 0.25 solar value [3l[. If the dense cloud we consid- 
ered has such a low metal abundance, the cosmic-rays 
from the cloud does not increase the observed B/C ra- 
tio. Moreover our scenario requires that the dense cloud 
has not been destroyed by stellar winds before the ex- 
plosion of the stars. It is known that stellar winds are 
weaker for stars with lower metal abundances, and that 
the metal abundances of stars are correlated with those 
of the host cloud. Thus the cloud in which electrons and 
positrons are effectively generated may tend to have low 
metal abundances. In summary our scenario cannot be 
rejected only by the present observational data of the 
B/C ratio. 
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